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We investigate stabilizer codes with carrier qudits of equal dimension D, an arbitrary integer 
greater than 1. We prove that there is a direct relation between the dimension of a qudit stabilizer 
code and the size of its corresponding stabilizer, and this implies that the code and its stabilizer are 
dual to each other. We also show that any qudit stabilizer can be put in a standard, or canonical, 
form using a series of Clifford gates, and we provide an explicit efficient algorithm for doing this. 
Our work generalizes known results that were valid only for prime dimensional systems and may be 
useful in constructing efficient encoding/decoding quantum circuits for qudit stabilizer codes and 
better qudit quantum error correcting codes. 
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PACS numbers: 03.67.Mn, 03.67.Pp 



I. INTRODUCTION 

Quantum error correction is an important part of 
various schemes for quantum computation and quan- 
tum communication, and hence quantum error correcting 
codes, first introduced about a decade ago [lHl] have re- 
ceived a great deal of attention. For a detailed discussion 
see Ch. 10 of Most of the early work dealt with codes 
for qubits, with a Hilbert space of dimension D = 2, but 
qudit codes with D > 2 have also been studied [5Hl5j. 
They are of intrinsic interest and could turn out to be of 
some practical value. 

The stabilizer formalism introduced by Gottesman in 
16] for D — 2 (qubits) provides a compact and pow- 
erful way of generating quantum error correcting codes 
and extends the notion of linear classical error correcting 
codes (l7j to the quantum domain. The stabilizer for- 
malism has been generalized to cases where D is prime 
or a prime power, see e.g. [f| [H, [H, EH- For composite 
D things are more complicated and there is no imme- 
diate and natural way of generalizing the notions. Our 
approach is to use generalized Pauli operators and stabi- 
lizers defined in the same way as in the prime case, see 
e -g- [HI EH]- This has the virtue that many (although 
not all) results that are valid in the prime dimensional 
case can be extended without too much difficulty to the 
more general composite case. 

An important problem in the theory of stabilizer codes 
is what is their structure. Is there any "canonical" way 
of representing an arbitrary stabilizer code? If yes, can 
one use this fact for implementing various quantum error- 
correcting tasks? For prime D it turns out that there is 
such a standard form, see e.g. Ch. 10.5.7 of [H, and this 
allows for a better understanding of the error-correcting 
capabilities of the stabilizer code and also provides an 
efficient way of constructing encoding/decoding circuits 
for such stabilizer codes. For composite D we are not 
aware of any such standard form (except for the case of 



stabilizer codes over prime-power finite fields [12|). and 
the proof that such a form exists is one of the main results 
of the current article. 

The reminder of the paper is organized as follows. 
Sec. |TT] contains definitions of the generalized Pauli group 
and some quantum gates used later in the paper. It 
also defines rigorously qudit stabilizers and their corre- 
sponding stabilized subspaces (or codes), together with 
an alternative algebraic notation that we employ later. 
Sec. Mil contains our main results: a "size" theorem that 
relates the size of the stabilizer group to the dimension of 
its stabilized subspace, followed by a "structure" theorem 
that shows that any qudit stabilizer can be brought to 
a standard form through a series of elementary quantum 
gates. Finally, Sec. II VI contains a summary, conclusions, 
and some open questions. 



II. PRELIMINARY REMARKS AND 
DEFINITIONS 

A. The generalized Pauli group on n qudits 

We generalize Pauli operators to higher dimensional 
systems of arbitrary dimension D in the following way. 
The X and Z operators acting on a single qudit are de- 
fined as 
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and satisfy 



X D = Z D = I. 



XZ = uZX, 



(1) 



(2) 



where the addition of integers is modulo D, as will be 
assumed from now on. For a collection of n qudits we use 
subscripts to identify the corresponding Pauli operators: 
thus Zi and Xi operate on the space of qudit i. The 
Hilbert space of a single qudit is denoted by H, and the 
Hilbert space of n qudits by H n , respectively. Operators 
of the form 
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will be referred to as Pauli products, where A is an integer 
in Z^j and x and z are n-tuples in Z^ , the additive group 
of n-tuple integers mod D. For a fixed n the collection 
of all possible Pauli products ([3]) form a group under 
operator multiplication, the Pauli group V n . If p is a 
Pauli product, then p D — I is the identity operator on 
"H n , and hence the order of any element of V n is either 
D or else an integer that divides D. While V n is not 
Abelian, it has the property that two elements commute 
up to a phase: 



V1V2 = w Xl2 p 2 pi, 



(4) 



with A 12 an integer in Z^j that depends on pi and p 2 - 



B. Generalization of qubit quantum gates to 
higher dimensions 



Pauli operator 


S q 


F 


Z 


z q 


X 


X 




Z D-1 



TABLE I. The conjugation of Pauli operators by one-qudit 
gates F and S q (q is the multiplicative inverse of q mod D). 



Pauli product 


CNOT ab 


SWAP ab 


CP a j, 


la <g> Z b 


Z a ® Z b 


Z a ®I b 


la ® Z b 


Z a ® I b 


Z a ®I b 


la ® Z b 


Z a ®h 


I a ®X b 


Ia®X b 


X a ®I b 


ZP~ x ®X b 


X a ®I b 


x a ® xP- 1 


la ®X b 


Xa ® Z°~ X 



TABLE II. The conjugation of Pauli products on qudits a and 
b by two-qudit gates CNOT, SWAP and CP. For the CNOT 
gate, the first qudit a is the control and the second qudit b 
the target. 



In this subsection we define some one and two qudit 
gates generalizing various qubit gates. The qudit gener- 
alization of the Hadamard gate is the Fourier gate 
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F:=-^L$>' fe |j) (k\ 
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For an invertible integer q £ (i.e. integer for which 
there exists q € Z^ such that qq = 1 mod D), we define 
a multiplicative gate 



D-l 

S q :=£>> (jq\, 



(6) 



where qj means multiplication mod D. The requirement 
that q be invertible ensures that S 9 is unitary; for a qubit 
S q is just the identity. 

For two distinct qudits a and b we define the CNOT 
gate as 

D-l D-l 
CNOT ab := li) V\a® X i = E b') 01a®l*) ( fc + A . 

(7) 

the obvious generalization of the qubit Controlled-NOT, 
where a labels the control qudit and b labels the target 
qudit. Next the SWAP gate is defined as 



with 

(CNOT ba )t = (CNOT^- 1 = (7 Q ®F2)CNOT ba (/ a ®F2) 

(10) 

Finally we define the generalized Controlled-phase or CP 
gate as 

D-l D-l 
j=0 j,fc=0 

(11) 

The CP and CNOT gates are related by a local Fourier 
gate, similar to the qubit case 

CNOT ab = (J a ® F b )CP a6 (/ Q ® F 6 )t, (12) 

since F maps Z into A under conjugation (see Table [XJ) . 

The gates F, S„ SWAP, CNOT and CP are unitary 
operators that map Pauli operators to Pauli operators 
under conjugation, as can be seen from Tables Ulandllll 
They are elements of the so called Clifford group on n 
qudits [IFJ, [2l|, the group of n-qudit unitary operators 
that leaves V n invariant under conjugation, i.e. if O is 
a Clifford operator, then \/p e V n , OpO^ G V n . From 
Tables Q] and [TT] one can easily deduce the result of con- 
jugation by F, S g , SWAP, CNOT and CP on any Pauli 
product. 



D-l 



SWAP Qb := I*) 01„ ® li) (k\t 

j,k=0 



(8) 



It is easy to check that SWAP gate is hermitian and does 
indeed swap qudits a and b. Unlike the qubit case, the 
qudit SWAP gate is not a product of three CNOT gates, 
but can be expressed in terms of CNOT gates and Fourier 
gates as 

SWAP ab = CNOT ab (CNOT bQ )t C NOT Qb (F2 ® I b ), (9) 



C. Qudit stabilizer codes 

Relative to this group we define a stabilizer code C to 
be a K > 1-dimensional subspace of the Hilbert space 
satisfying three conditions: 

CI: There is a subgroup S of V n such that for every s 
in S and every \ip) in C 

<#) = m (13) 
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C2: The subgroup S is maximal in the sense that every 
s in V n for which (| 13[) is satisfied for all \ip) E C 
belongs to S. 

C3: The coding space C is maximal in the sense that any 
ket \4>) that satisfies (fl3|) for every s € S lies in C. 

If these conditions are fulfilled we call S the stabilizer 
of the code C. That it is Abelian follows from the commu- 
tation relation (j4j , since for K > there is some nonzero 
\tp) satisfying ([TBI . 

Note that one can always find a subgroup S of 
satisfying CI and C2 for any subspace C of the Hilbert 
space, but it might consist of nothing but the identity. 
Thus it is condition C3 that distinguishes stabilizer codes 
from nonadditive codes. A stabilizer code is uniquely 
determined by S as well as by C, since iS determines C 
through C3, so in a sense the code and its stabilizer are 
dual to each other. 



D. Stabilizer generators and equivalent algebraic 
descriptions of qudit stabilizer codes 

Any stabilizer group can be compactly described using 
a set of group generators. A generator corresponds to a 
specific Pauli product and can be completely specified, 
see ©, by a phase A and two n-tuples in Z^, x and z. 
A collection of k generators can therefore be represented 
by a fc-component phase vector over Z/j (that contains 
all k phases) and a k x 2n parity- check matrix over Z^> 
with rows corresponding to the stabilizer generators. For 
example, the stabilizer 

S=(u, 2 XfZl,Xi) (14) 

;ponds to the phase vector (2, 0) and parity-check 



corres 
matrix 



S = 













: 


2 








(15) 



The angular brackets in (|14[) means "group generated 
by" , i.e. the group obtained by all possible products of 
the group generators. We call the left k x n block of the 
parity-check matrix the A-block, and the right k x n the 
Z-block, since they describe the X and Z parts of the 
stabilizer generators, respectively. 

Note that if D is a prime number, any stabilizer group 
can be described using no more than n generators. How- 
ever, in composite dimensions one can have more that n 
generators but no more than 2n. For example, in D = 4, 
the n — 1 qudit stabilizer S = (X 2 , Z 2 } is generated by 
2 (and not 1) elements and specify the stabilizer state 
(|0> + |2))/V2. There is no way of representing this state 
using only 1 generator; Z 2 by itself stabilizes both |0) 
and 1 2), hence everything in their span, so the condition 
C3 is not satisfied, i.e. the coding space is not maximal. 
The same kind of analysis holds for X 2 . A more rigorous 



Gate 


X-part 


Z-part 


SWAP at 


Interchange columns 
a and b 


Interchange columns 
a + n and b + n 




Multiply column a by 
invertible integer g _1 


Multiply column a + n 
by invertible integer q 


(CNOT afc ) m 


Substract m times 
column a from column b 


Add m times column 
b + n to column a + n 



TABLE III. Conjugation by the above quantum gates corre- 
spond to elementary column operations on the X and Z parts 
of the parity-check matrix of a stabilizer code. For the CNOT 
gate, the first qudit a is the control and the second qudit b 
the target. The integer exponent m means CNOT applied m 
times (or, equivalently, the m-th power of the CNOT gate). 



analysis can be done using the Theorem Q] of Sec. IIII1 
which implies for this example that the size of the stabi- 
lizer group must be equal to 4, hence it must be generated 
by a single generator of order 4 or two generators each 
of order 2. By inspection it is easy to rule out the first 
case, so indeed one must use 2 generators. 

A conjugation of a stabilizer group by a Clifford op- 
eration will change the stabilizer group to an isomor- 
phic group. This will correspond to a column opera- 
tion on the parity-check matrix of the stabilizer, together 
with a transformation of the phase vector. On the other 
hand, the generator description of a stabilizer group is 
not unique: one can multiply a generator by another one 
and still get the same group. This kind of operation cor- 
responds to a row operation on the parity-check matrix, 
again keeping in mind that in general the phase vector 
will modify. From now on for the sake of simplicity we 
will ignore the phase vector, although in real applications 
one has to keep track of the phases. 

The following represent what we call elemen- 
tary row/column operations: a) interchanging of 
rows/columns, b) multiplication of a row/column by 
an invertible integer, c) addition of any multiple of a 
row/column to a different row/column. The column op- 
erations can be realized by conjugations of the stabilizer 
by the Clifford operations in Table IIII1 and the row op- 
erations just ensure that the stabilizer group remains the 
same, i.e. the new set of generators generate the same 
stabilizer group and not a smaller one. 



III. SIZE-STRUCTURE THEOREMS 

The following theorem generalizes to composite D a 
well-known result for prime D that relates the size of the 
stabilizer group to the dimension of its stabilizer sub- 
space. Although the composite D result of our next the- 
orem may have been known by the community (see e.g. 
the claim near the end of Sec. 3.6 of [HI), we have not 
yet seen a proof of it. 

Theorem 1 (Size). Let C be an n-qudit stabilizer code 



4 



with stabilizer S. Then 



K x\S\= D r 



where K is the dimension of C, \S\ is the size (or order) 
of the stabilizer group S and D is the dimension of the 
Hilbert space of one carrier audit. 



Proof. Define 



P 



151 2' 



(17) 



ses 



We will first show that P is the projector onto C. 
It follows at once that 

p = pt = 



where the equalities follow from the group property of iS, 
so P is an orthogonal projector. Let be an arbitrary 
vector that belongs to the stabilizer code C. Then s\ip) = 
Vs (see the condition CI (TIU]) that a stabilizer code 
must satisfy), which together with (fT?]) yields 



P\tf>) = |# 



the 1-qudit stabilizer generated by S = (Z 2 ) in D = 4. 
It is obvious that S stabilizes a K = 2-dimensional code 
(16) C=span{|0), |2)}, and the size of the stabilizer is |«S| = 2. 

Whenever D = 2 (qubits) it was shown in [l6| (see also 
Ch. 10.5.7 of [4J for a detailed discussion) that any stabi- 
lizer code can be put into a "standard form" or "canoni- 
cal form" , and this is very useful for constructing encod- 
ing/decoding quantum circuits for stabilizer codes. This 
result can be generalized at once to prime D. However, 
for composite dimensions, it is not so obvious how to do 
the generalization, and the main technical difficulty is 
that Zj3 is a ring (and not a field) and therefore some 
integers are not invertible. However, in the following 
Theorem we show that one can still apply a technique 
similar to a Gaussian elimination over and put any 

(18) composite D stabilizer code into a standard form similar 
to the one of prime D case. 

Theorem 2 (Standard form). Let C be a K dimen- 
sional n-qudit stabilizer code with stabilizer S generated 
by k In generators and with corresponding parity-check 
matrix S of size k x 2n. Then S is isomorphic through a 
conjugation by a Clifford operation to another stabilizer 

(19) S' , with parity check matrix S' in standard form 



Therefore the subspace W onto which P projects includes 
C, C C W. Let us now choose an arbitrary \(f>) € W. 
Then 



PL 



151 2' 



(20) 



ses 



Multiply ([2"0")l on the left by some arbitrary t £ S and use 
the group property of S to get 

t|*>=ii'i0> = ii|5>w = t^E s I^ 

' ' seS ' ' sS5 

= p\4>) = \4>). (21) 

Since t was arbitrary we arrived at the conclusion that 

t|0) = |0),VieS, (22) 

which proves that \<f>) belongs to the stabilizer subspace 
C, and this implies that W C C. Hence P projects strictly 
onto C (and not on some larger subspace that includes 
C). Its trace is just the dimension K of C, 



Tr(P) = K= ifl 1 



(23) 



where we have used that fact that all Pauli products in 
(fl"7|) are traceless except the identity (that must belong to 
the sum, since S is a group), of trace D n . This concludes 
the proof. □ 

When I? is a composite integer the dimension of the 
stabilizer code does not have to be a power of D any 
more (as was the case in the prime D case), but can be 
any divisor of D n . As an illustrative example, consider 



S' 





(24) 



where the dimensions of the block matrices are indicated 
by curly brackets. 

Here M = diag(m\ , . . . , m r ), with 1 ^ r ^ n, is a diag- 
onal matrix with all rrij 7^ divisors of D. The matrices 
Z\ and Z2 satisfy 



Z 1 M= MZ[ 
Z 2 M= 



mod D, 
mod D, 



(25) 
(26) 



where T in the exponent denotes the transpose, and the 
matrix Z^ is a diagonal rectangular matrix, with diagonal 
elements divisors of D. The notation denotes the zero- 
block matrix. 

Proof. The key ingredient of the proof is the Smith nor- 
mal form: through a sequence of elementary row/column 
operations mod D (see the discussion at the end of 
Sec. IIII1 a matrix M over can be converted to the 
Smith normal form [22|,[23[ (see also Sec. IV. B of [l5[ for 
an example) 



M' = V M W, 



(27) 



where V and W are invertible (in the mod D sense) 
square matrices, and M' is a diagonal rectangular ma- 
trix, with diagonal elements divisors of D. The matrix 
V represents the row operations and W the column op- 
erations. 

Note that in our case all necessary column operations 
can be realized by the corresponding gates in Table IIII1 
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and, more important, they do not mix the X and Z parts 
of the parity-check matrix. Therefore, without being con- 
cerned with what happens to the Z part of the parity- 
check matrix S, we can put its X part in the Smith nor- 
mal form (again we stress that this can be done because 
the Z part to not interfere with the X part), and ar- 
rive at a parity-check matrix of the form (|24l) . Next by 
another series of Clifford gates acting only on the last 
(n — r) qudits one can further put the Z4 matrix in its 
Smith normal form, without modifying the X part of 
the parity-check matrix (which is already in Smith nor- 
mal form), since there are only zeros on the last n — r 
columns of the X part; note that the row operations are 
done on the last k — r rows, and again do not modify the 
X part of the parity-check matrix. Since the elementary 
row operations do not change the stabilizer group and 
the elementary column operations correspond to Clifford 
gates, see Table Hill our whole transformation from S to 
S' is a conjugation by a Clifford operation. 

Finally note that conjugation by Clifford operations 
do not change the commutation relations. It is easy to 
deduce from © that two Pauli products described by 
(x|z) and (x'|z') commute if and only if 

x z' = z • x' mod D, (28) 

where the dot represents the usual inner product of two 
vectors in Zb, e.g. the sum of the products of individual 
components. Using (|28[) we observe at once that the final 
set of generators commute if and only if (|2~5|) and (|26l) 
hold, and this concludes the proof. □ 

It is proved in [23[ that a M x N matrix can be re- 
duced to the Smith form in only (D(M N) operations 
from TL-q , where 6 is the exponent for matrix multiplica- 
tion over the ring Z^, i.e. two M x M matrices over Z^ 
can be multiplied in 0(M ) operations from Z^. Using 
standard matrix multiplication # = 3, but better algo- 
rithms [24j allow for 9 = 2.38. This ensures that our 
procedure outlined above is computationally efficient. 

Whenever D is prime its only non-zero divisor is 1, 
hence M is just the identity matrix and Z4 has only O's 
and l's on the diagonal. From (|25l) and (l2l)|) . Zi must 
be symmetric and Z2 must be the zero matrix, hence our 
standard form reduces to the one known for prime Z?'s 
16]. 



Also for prime D it was proven in [12| that any stabi- 
lizer group S is Clifford equivalent to another stabilizer 
S' generated only by Z% S' = {Z x , Z 2 , ■ ■ . , Z k ) 0. This 
result is not true for composite -D's, and one easy to see 
counterexample is the 1-qudit stabilizer S — (X 2 , Z 2 ) in 
D = 4, mentioned before in Sec. Ill Dl 
IV. CONCLUSIONS AND OPEN QUESTIONS 

We studied stabilizer codes with carrier qudits of com- 
posite dimension D. We proved a size theorem that re- 
lates the size of the stabilizer group to the dimension of 
its stabilized code. Furthermore, we have shown that 
any stabilizer code can be put in a standard (canonical) 
form through a series of Clifford gates, and we provided a 
constructive algorithm. Our result generalizes what was 
known in the prime D case and may be useful in con- 
structing efficient encoding/decodin g q uantum circuits, 
following the procedures outlined in [16j. 

Our approach was based on the generalized Pauli group 
introduced by ([U and ((2|). However, for composite di- 
mensions, this is not the only way of introducing Pauli 
operators. An alternative way is to split the dimension in 
its prime-power factors which will then induce a natural 
splitting of the carrier qudits in subsystems of prime- 
power dimensions. In each of these subsystems then one 
can define Pauli operators using finite fields (any finite 
field is isomorphic to a prime-power canonical represen- 
tation), as done e.g. in 12]. Although this is the scope 
of future work, we think it may be useful since in a sense 
"decouples" the stabilizer into prime-power subsystems, 
and the latter can be put into standard forms as done 
in [l2j | . One can then use previously known results for 
stabilizers over finite fields to study various properties of 
composite D stabilizers, and this may help building more 
efficient quantum error correcting codes. 

Finally one may ask if there exist alternative standard 
forms of qudit stabilizer codes, perhaps more useful that 
the one presented here. We do not know if such forms 
exist, and searching for them may be worthwhile. 
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